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ABSTRACT 


The coefficient of the PIGA second-order nonlinearity term is 
developed in terms of the misalignment between the PIG and PIGA axes 
and the misalignment of the PIG axes with respect to the principal axes 
of inertia. The results indicate that to further reducé the nonlinearity 
term, ao, equalization of the moments of inertia about SRA and IA is not 
profitable. Emphasis should be placed on dynamically balancing each of 
the rotating components, with a view towards reducing center of gravity 
shifts which contribute to the misalignment angle between the PIG IA and 


the principal axis. 
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STUDY OF THE SECOND-ORDER TERM 
IN THE MOD J 16 PIGA 


1, INTRODUCTION 

The intent of this study is to define the contributions to the second-order term in the 
Mod J 16 PIGA. Several reports, as listed in the references, have dealt with torques due to float- 
wheel misalignments, mass unbalances, and flotation, but in view of the rigid specifications im- 
posed on the Mod J 16 PIGA, a definitive study of second-order effects was undertaken, There- 
fore, the complete second-order term of a nonideal PIGA operating under a constant acceleration 


will be defined. 


2, DEFINITION OF THE SECOND-ORDER TERM 


The acceleration indicated by a PIGA output expressed in a power series is: 


where: 
Oe indicated PIGA acceleration output 
Oe is acceleration input 
Oo» 81> Bor a, 4° * = coefficients 


The coefficient a5 is known as the second-order term. The torque summation equation about 


the PIG output axis under steady-state conditions is: 


(2) 


where: 


| H Wa =~ gyroscopic torque about OA 
H = angular momentum of gyro wheel 
Win = angular velocity of PIG about IA 


AME TOONS ARON NMA eee 


mé a, = pendulous torque about OA 
m = _ pendulous mass 


f = pendulous arm 


a. = component of acceleration along PIGA IA 


M, = inaccuracy torque 


In this study, eenaideraion is given to the inaccuracy torques dependent on the square of the 
applied acceleration only. Therefore, since the angular velocity of the PIG Eeat IA is propor - 
tional to the applied acceleration, let: 


eS 2 ; 
M, = AWA (3) 
where: 
X% = constant of proportionality 
Substituting Eq. (3) into Eq. (2): 
2 & : 
AWra - HWy, + még ain 0 (4) 


Solving the quadratic equation (4) for Wy: 


WwW Get oily (5) 


IA 2X 2X 


The only solution that describes the physical rotation of the PIG about IA is: 


(6) 


The other solution results in an alternating series expansion for Wra: Factoring Eq. (6): 


(7) 


Using an expansion series: 


DS ye ee ee et (8) 


and letting x = (4mé4 a,/H) (A/T), Eq. (7) becomes: 


2 \2 A 
= ee (_) ta (10) 


Since the parameters m, £, and H are defined by design, the parameter A is the only undefined 


parameter required to define the second-order term, 


3. DEVELOPMENT OF INACCURACY TORQUES AS A RESULT OF PIG TO PIGA MISALIGNMENT 


Hook ' 
Consider the orthogonal set of axes 0x y z shown in Fig. 1 where the Ox axis is defined 
as the PIGA input axis and where the angular velocity of the float about the PIGA IA is W. The 
total torque about the point 0 as a result of the angular velocity about the PIGA IA is: 


M, = as R, x F, (11) 
where: 
M, = torque about point 0 
F; = centrifugal force acting on the mass m; 


distance of mass m,; from point 0 


- etneontinn: 


Y= TOTAL MISALIGNMENT ANGLE BETWEEN PIG IA 
AND PIGA JA 


Y, * COMPONENT OF PIG TO PIGA IA MISALIGNMENT ANGLE 
ABOUT SRA IN THE IA-OA PLANE 


Yo = COMPONENT OF PIG TO PIGA IA MISALIGNMENT ANGLE 
ABOUT OA IN THE IA-SRA PLANE 


Fig. 1 Components of the PIG te PIGA Misalignment Angle. 


_ ‘The centrifugal force, Fi, is: 
F, =m; r; W (12) 


! 
where r; is the distance of the mass m, from the 0x axis and the summation is taken for all 


particles of mass m, that make up the Hoa structure. The component of the total torque about 
the oy’ axis isi 


i] 1 
' ' s ' 
M, a, (2; Pix yr ak a) 
where: 
! 
ey = component of force EF; acting along the Ox axis = 0 
: ‘ 
' 
Bae’ = component of force F, acting along the 0z axis 
Pee ei 


' 
component of vector R; along the 0x axis 


3 ' 
The component of the force acting along the 0z axis is: 


Ri) 8 F, cos (30-6,) (14) 
where: 
H 
9; = angular location of Fy from the Oy axis. 
" Substituting Eq. (12) into Eq. (14): 
Fi! ys mi," e, Wo ein €, (15) 
1Z a i 1 


Substituting Eq. (15) into Eq. (13) and,since F;,' equals zero, the torque about the Oy' axis 
becomes: 


M,,’ = mua > x, m; Yr; sin 0; (16) 


| 
H 
i 
| 
i 
/ 
i 


pd ivies 


rer bene dmtndaee. 6 


atesnias 


i 
Since Raa sin 9;, Eq. (16) becomes: 
x 2 ' ' : 
M = -W du. x; 2; mM, (17) 


. ! f 
Theretere, the dyAamis uaApalance fAFAHB AbBUt the OY axis is: 


y xz (18) 


where: 


fat i 
I,',' = float product of inertia in the x z plane = ~- ‘S m, x. Zz. 
ies a 

The torque about the float OA axis (Oy axis) is: 


EON eats My! S28 Arnie OA) ue) 


where: 
Moa = torque about float OA, which is also defined 
as the Oy axis 


i 
A = misalignment angle between the Oy and Oy axes 
(OA misalignment angle) 


Neglecting the float OA misalignment angle: 
Moan = M SV Gee Wee elt (20) 


The float OA misalignment angle can be neglected in Eq. (19) because an OA misalignment angle 
that is © 20 arc-sec, for example, results in an error of less than 0,05 ppm. To express the 
float product of inertia ae in terms of the float coordinate frame IA, OA, SRA earani ra through 
the PIG to PIGA misalignment angle y (see Fig. 1), consider a matrix composed of two 
consecutive rotations: 


1 
clockwise rotation about the 0z axis 


(1) ¥, 


' 
(2) Y9 clockwise rotation about the Oy axis 


The eeerdinates ef M; expressed in the transformed frame are: 


x cos Yo 0 sin Yo COs ¥, sin 14 0 x 
' 

y = 0 ] 0 -sin ¥,; cosy, 0 y (21) 
' 

", -sin Yy 0 cos V9 0 0 1 Zz 

where x, y, z are the distances along the float IA, OA, and SRA axes, respectively, 
. 

' 

x X COS Yo Cos v; + y sin 1 COS Yo + z sin Yo 
' 

y | =| -x sin ¥, + ycos 11 (22) 
' 

Zz -x sin Yo cos >, he ay sin Y; sin Yo + Z cos Yo 


Using small angle approximations where the cos of an angle equals one and the sin of an angle 
equals the angle, Eq. (22) becomes: 


x 2 eS Vay ov: Yo% 
i 

a. = “7, es iy (23) 
' 

: at eis eS 


From Eq. (23), the product of x 2. is: 


dara a 3g 2 
ee ae aaa (27172) eet (: : 7) x2 + (rg * 7) yz + 92” : ¥ oy" (24) 


Neglecting products of angles, Eq. (24) may be expressed as: 


° 


ti 2 
XZ = Yo {c + y*) 2 (x? + ”)]. Y, yz + xz (25) 


3 


enema 


[a | 
The product of inertia in the x z plane is: 
' ! 
ee Ears ae 
ar > ses ae io (26) 
Substituting Eq, (25) into Eq. (26): 


2 7, 2 
T'2' = Ye dom (x? + 4) - Te dp mte” +y") - y S mye Domxz, M27) 


Tiges® g(t, * Ee) ty I. T Te ? (28) 
where: . 
we = moment of inertia about the Ox axis 
L. = moment of inertia about the 0z axis 
le = product of inertia in the yz plane 
le = product of inertia in the xz plane 


Based on Eq. (20), the torque about the float output axis is: 


2 


Moa = ee (29) 


Substituting Eq. (28) into Eq. (29), and neglecting Y 1 with respect to the other terms, the por- 


tion of the error torque due to mass unbalance and due to the misalignment between the PIG and 


PIGA input axes is: 
M © Le Sep our Fee (30) 
OA YQ ZZ xx xZ 


4, DEVELOPMENT OF PIGA NONLINEARITY COEFFICIENT AS A RESULT OF PIG TO PIGA 
MISALIGNMENT AND ALIGNMENT TO FLOAT PRINCIPAL AXES 


Consider the three coordinate systems shown in Fig. 2. The coordinate system xyz 
is an orthogonal set of axes with the Ox: axis defined as the PIGA IA. The coordinate system 
xyz is an orthogonal set of axes that defines the float IA, OA, and SRA, respectively. The 
coordinate system Xp Le zy is an orthogonal set of axes that defines the principal axes of inertia 
of the float. 


A rigid body rotating about a principal axis of inertia and not subject to an external torque 


tr 


2p 


x = 
f°) y 
PIG 
p OA 
it~ Y} % 
Bo By x! 
PIGA 
lA 


Fig. 2 Coordinate Systems Representing the PIGA Axes, the PIG Axes, 


and the Principal Axes of inertia, 


will rotate in equilibrium. In addition, 
when rotation is about the principal axis 
16 PIGA when: 


Stable equilibrium of a Symmetrical body occurs only 
of greatest moment of inertia, This is true of the 


fa? leg a "Toa (31) 
where: 
a = moment of inertia about float IA 
Iona = moment of inertia about float SRA 
Toa * moment of inertia about float OA 


‘ 
a centrifugal couple, proportional to the 
ng minember, 


When the axis of rotation is not about a principal axis, 


square of the angular velocity, is exerted on the rotati 


The product of inertia iS! expressed in terms of the xyz set of orthogonal axes is given 


in Eq. (28), In like manner, the product of inertia ee expressed in terms of the orthogonal axes 


representing the principal axes of inertia, 


Xp¥ pp’ is: 
I, 5 B, (1, i 1, } - ae (32) 
p p Pp 
where: 
By = rotation angle about the Oy axis 
I = moment of inertia about 0z 
P = 
a. = moment of inertia about 0x 7 
p 
XZ = product of inertia in the %>%p plane 
Since by definition of principal axes, all products of inertia vanish: 
I I aan) (33) 
xX "2 
P Pp 
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Substituting Eq. (32) into Eq. (28), and assuming that (1, = ,) 3 (4 ar ), Eq. (28) 
becomes: P P 


Ie " i 7 Ses . Bs ) (34) 


3 (sna : ta) (2 bl By] Pe 


- where: 


Yo = component of PIG IA to PIGA IA piceeUrnment 
angle in the IA-SRA plane 


B = component of PIG IA to principal axis misalignment 
2 
angle in the IA-SRA plane. 


Substituting Eq. (35) into Eq. (29), the torque about the float or PIG OA is: 


% . 2 
Moa (sa ia) (1 : Pa) Wia 86) 
Substituting Moa = H Wia and we, = (SF)? a? n into Eq. (36): 
HW,, = (SF)? |(1 as + Bs || 22 (37) 
IA SRA TAY \¥2 2/| *in 
Referring to Eq, (10), the term 2 is: 
A = (Tee . ha) (7 + B, | (38) 
and the second-order term is: 
Bg (sr)? A 
' “ ee (Ispa lia \(v, + B, ) rad/sec 
2 H 2)2 
cm/sec 


In terms of the more familar Fy coefficient: 
1 


I eae) Yo + 
F, = (g) (sr) USRA ws) (% x 10° ppm/g (40) 
| H 
whe Fei 
g = 980,39215 cm/sec” 
SF = 0;00102 (rad/sec)/(em/sec?) 


5. METHOD FOR ESTIMATING COMPONENT Ol* MISALIGNMENT ANGLE BETWEEN PIG IA 
AND PRINCIPAL AXIS OF INERTIA IN THE IA-SRA PLANE 


In the expression of the second-order term, Eq. (39), the terms SF, Iona? Tae and H 
are defined by design, and Yo: the component of the total input axis misalignment angle y in the 
IA-SRA plane, is adjusted to a minimum during assembly (see reference 5). A method of estimat- 
ing the remaining term in Eq. (39), Bo, is described herein. The sequence of steps required to 
estimate By is as follows: 


(a) Measure Ioa? lia: and Iona of the float, 


(b) Estimate Iopa-oa’ Ispa-rae and Tia-oa Of the float, 


(c) Form the float inertia matrix Ip 


(d) Subject the float inertia matrix toa principal axis transform, 


(e) Obtain the direction cosines of the principal axes of inertia 


Xiunyez lati 
pp p re ative to the float axes xyz. 


(f) Calculate By = 1/28, 


(1) The moments of inertia about the float axes, Ina: Ta? and Iona can be measured 


to an accuracy of + 0.04%, using a torsional pendulum (see reference 6), 


(2) The products of inertia IcRA-OA? Top A -IA? and LA-OA can be obtained only = 
estimating the minimum and maximum values. The values used were 0.001 and 0,1 gm-cm", 


(3) The float inertia matrix is formed as: 


Toa Toa -1A TOA-SRA 
nes Joa -1A Ta Ign a-1A (41) 
Toa-SRA TsrA-1A ° Igpa 


(Since Ls is composed of coefficients of inertia, I, is symmetric.) 


12 


(4) The principal axis transform consists of diagonalizing I 
of inertia about the principal axes. To diagonalize I 
from the determinant 


Fr? giving the three moments 
pe Solve the cubic equation in I resulting 


Togs! Toa -1A Toa-SRA 
Toa-1A Taw! TsrA-1A | = 9 (42) 
ToOA-SRA Ispa-tA— [gral 


The three roots are the three principal moments of inertia I L I 9 and I 3 
p p3" 
(5) The direction cosines are obtained from the matrix transformation that transfers the 


xyz axes into the alignment with XD pp axes. This is obtained by solving the following simultaneous 


equations for the ratios y/x and z/x, ; 


(Fos - Ipi) * Toa-a¥ + Ipa-spa2 = 0 4s) 
lance © (4a - Tp) ¥ dae eae Pe te an (44) 


The result is normalized by Setting x = 1,0 and the angle B is calculated from: 


cosB = ——"—________. (45) 


| (6) The angle By is calculated as: 
me ET: (46) 


6. CONCLUSIONS : 


As indicated by Eq. (40), only the components of the misalignment angles y and 8 in the 
SRA- IA plane contribute to Fj,. In all the calculations, the assumptions were made that B ae 
1/2 B and Ct Wee 1/2 y. A computer program was written for the sequence given in paragraph 5 
to determine £B 2 for various PIGA units with products of inertias varying from 0.001 to 0.1 


gm-cm?, The results of this program and the calculations of the Fy 
1 


term are given in Table 1, 
Fig. 3, and Fig. 4, 
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Equation (40) seems to indicate that mere equalization of the moments of inertia about 


SRA and IA is sufficient to control the Fy term, This assumes that Bo remains constant as Ion A 


1 
and la are changed to reduce their difference, This assumption is incorrect because Bo increases 


rapidly as the difference between Ton A and La decreases (see Fig, 3), As shown in Fig, 3, a de- 

crease in the difference of Iona and La from 2,015 to 0,154 gm-cm? results in a decrease in the 

hy term 6f dii average of only 0:075 ppm/g: For Ioa-id 29d IshA-GA less than 0; 01 em -en? 

(which is probably more correct), the decrease in the Fy term is much smaller than 0,075 ppm/a. 
] 


The conclusion is reached that the main effort in reducing the F; term should be di- 
1 
rected toward dynamic balancing of each rotating component to reduce the shifts in the center 


of gravity that contribute to the angle Bo, (see reference 6). Further reduction in the inertia 


difference between Ign aq aid Ip, is not profitable due to the relatively little effect on the F, term, 
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